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ABSTRACT 

Rhizopus oryzae is a filamentous fungus, belonging to the order Mucorales. It can ferment a wide 
range of carbohydrates hydrolyzed from lignocellulosic materials and even cellobiose to produce 
ethanol. However, R. oryzae also produces lactic acid as a major metabolite, which reduces the 
yield of ethanol. In this study, we show that significant reduction of lactic acid production could be 
achieved by short (25nt) synthetic siRNAs targeting the IdhA gene. The average yield of lactic acid 
production by R. oryzae during the batch fermentation process, where glucose had been used as a 
sole carbon source, diminished from 0.07gm/gm in wild type to O.Olgm/gm in silenced samples. In 
contrast, the average yield of ethanol production increased from 0.39gm/gm in wild type to 
0.45gm/gm in silenced samples. These results show 85.7% (gm/gm) reduction in lactic acid 
production as compared with the wild type R. oryzae, while an increase of 15.4% (gm/gm) in 
ethanol yield. 

KEYWORDS: siRNA, delivery, Rhizopus oryzae, IdhA gene, lactate dehydrogenase, lactic acid, ethanol 
production 



INTRODUCTION 

Lignocellulosic biomass is considered as the prime 
component for use as carbohydrate feedstock for industrial 
fermentation of chemicals. Currently, ethanol is the largest 
source of biofuel in the global market, which is produced 
either from sugar-based materials, such as sugarcane, or 
from starch-based materials, such as corn. However, there 
is an intensive global research effort to develop the 
process of ethanol production from lignocellulosic 
materials. This process involves pre-treatment of lingo- 
celluloses, acid or enzymatic hydrolysis to extract sugars 
from cellulose and hemicelluloses, followed by their 
fermentation into ethanol. Enzymatic processes require 



two classes of enzymes, cellulase and P-glucosidase, as 
well as the microorganisms to produce ethanol: 

fj-galactosidase 

cellobiose ► glucose 

* \ 

cellulase microorganisms 
/ V 

cellulose ethanol 



Since optimum temperatures for enzymatic hydrolysis 
and fermenting organisms are generally different (45°C 
vs 30°C, respectively), the process is sometimes carried 
out by separate hydrolysis and fermentation (SHF). 
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However, SHF suffers from the end-product (glucose) 
inhibition; therefore, another approach known as 
simultaneous saccharification and fermentation (SSF) 
carried out at a compromising temperature (e.g., 37oC) is 
preferred. 

There are several microorganisms that can utilize 
lignocellulosic sugars and ferment them into numerous 
valuable products. R. oryzae is one such microbe; it is a 
saprophytic fdamentous fungus and can up take a wide 
range of sugars such as, mannose, xylose, glucose, and 
galactose (Edebo, 2000), and the cellulose unit monomer 
cellobiose (Karimi et al, 2006). It has been largely used 
as a host microorganism to produce lactic acid (Skory, 
2004a), but it can also produce a variety of other valuable 
metabolites, such as, gallic acid (Misro et al, 1997), 
lipase (Salah et al, 1994), protease (Tunga et al, 1999), 
cellulolytic enzymes (Amadioha, 1993), and ethanol 
(Abedinifar et al, 2009). It is a good alternative for 
ethanol production because of its tolerance to inhibitors 
in lignocellulose acid hydrolyzates (Karimi et al, 2006), 
its valuable material contents in the biomass, and its 
ability to grow at higher temperatures lowering the risk 
of contamination (Abedinifar et al, 2009). Thus, R. 
oryzae can utilize cellobiose quite well in the SSF 
processes: 

cellulase R. oryzae 

cellulose ► cellobiose ► ethanol 

However, ethanol yield by R. oryzae fermentation is 
relatively low because the fungus converts the sugar into 
lactic acid as a major by-product (Abedinifar et al, 2009). 
Therefore, it is conceivable to silence genes involved in 
lactic acid production to increase ethanol yield. 

In RNA interference (RNAi), dsRNA facilitates 
degradation of the homologous mRNA, thereby 
diminishing or abolishing gene expression (Fire et al, 
1998). However, despite numerous reports on RNA 
silencing in a variety of organisms, only a few species 
among fdamentous fungi, such as Neurospora crassa and 
Cryptococcus neoformans have been shown to have 
RNAi machinery (Cogoni and Macino, 1999). 
Furthermore, RNA silencing has also been shown in 
Magnaporthe oryzae, the rice blast fungus by Naoki 
Kadotani and colleagues, using enhanced green 
fluorescence protein (eGFP) as a model (Kadotani et al, 
2003). Skory (Skory, 2004b) tried several types of gene 
expression systems for R. oryzae to produce altered 
transcript that can form dsRNA of the IdhA gene. 
However, he was not able to achieve any significant 
reduction in lactic acid production for isolates containing 
short (20-25nt) synthetic IdhA RNAi in the expression 
plasmids. 

The aim of this study was to address two main questions: 
(i) Can direct delivery of siRNAs results in silencing of the 
IdhA gene in R. oryzae? (ii) Can silencing of IdhA gene 
favorably increase ethanol production in the fermentation 
process considering the biochemical pathway of ethanol 
production in R. oryzae? 



MATERIALS AND METHODS 

Fungal strain, growth media and spore solution 

R. oryzae CCUG 28958 (Culture Collection, University of 
Gothenburg, Sweden) was used in all experiments. The 
strain was maintained at 4°C on potato dextrose agar 
(PDA) slants (potato extract 4gm/l, glucose 20gm/l, agar 
15gm/l). For spore formation, a swab from the slant was 
taken and spread onto PDA agar plates. The plates were 
incubated at 30°C for 4 days. To prepare the inoculum, the 
agar plates containing sporulated fungi were washed with 
20ml sterile water, and 20ml/l spore solution used in 
fermentation. 

siRNA design 

We targeted the region of IdhA with the highest homology 
with IdhB to silence both genes simultaneously in order to 
reduce chances of any bypass activity of IdhB. We used 
BLASTN 2.2.24+ (Zhang et al, 2000) to determine the 
regions of IdhA gene with high homology to IdhB gene. 
SiRNAs were then designed using these sequences with 
BLOCK-iT™ RNAi Designer (Invitrogen, Carlsbad, CA). 
The following sequences were used: 

Sense: 5 ' GGAGGC AGGGC AGGC AG AU AUUGUU 
Antisense: 5 ' AACAAUAUCUGCCUGCCCUGCCUCC 

The target sequence was screened against the Broad 
Institute database (Ustilago Maydis Sequencing Project, 
Broad Institute of MIT and Harvard) to ensure that the 
sequences were specific for the target genes. 

siRNA oligonucleotides (Invitrogen, Carlsbad, CA) were 
stored at -20°C until used. Duplex siRNAs (80nmol) were 
resuspended in 800(il DEPC treated water (Invitrogen, 
Carlsbad, CA), and then diluted 1:5 in TE buffer (lOmM 
Tris-HCl, 20mM NaCl, ImM EDTA, pH 8). 

Fungal protoplast formation and transformation 

Protoplast preparation and transformation in fungi has 
many limitations due to the wide variety of cell 
composition in different fungal species. Every so often the 
lytic enzymes used in the process do not digest the cell 
properly (Nakayashiki, 2008). We used the CaCl 2 /PEG 
(polyethylene glycol) method to transform R. oryzae with 
siRNAs, which has been used for some fungal species, 
such as Magnaporthe oryzae, Neurospora, Aspergdlus, 
and Fusarium (Nakayashiki et al, 2005). This method was 
originally reported by Vollmer and colleagues (Vollmer 
and Yanofsky, 1986). Briefly, 150ml of semi-synthetic 
growth media was inoculated with 3ml of spore solution 
and incubated at 30°C and with shaking at 140rpm, for 3-4 
days. Then R. oryzae mycelia were harvested by filtration 
through a Buchner funnel containing filter paper 
(Whatman No. 1). The mycelia pad was transferred to a 
plastic tube and 10ml of digestion buffer (lOmM 
Na 2 HP0 4 , 1.2M MgS0 4 , lOmg/ml Lysozyme) (Sigma- 
Aldrich) was added per gram, wet weight of mycelia. The 
tube was incubated for 3hr at room temperature with 
gentle inversion on a horizontal shaker (New Brunswick 
Scientific, USA). The contents of the tube were passed 
through four layers of sterile gauze. The filtrate was 
centrifuged at 800xg for 5min (Megafuge 1 .OR, UK), and 
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the collected protoplasts were resuspended in 50ml STC 
buffer (1M Sorbitol, 50mM CaCl 2 , 50mM Tris-HCl, pH 8) 
(Nakayashiki, 2008). After another round of centrifugation 
as above the supernatant was discarded and the protoplasts 
were prepared. 

Fungal transformation was carried out by mixing 100(0.1 
protoplasts and 100m1 of 20(o.M siRNA solution. The 
mixture was incubated on ice for 15min. 1.2ml of PEG 
solution (60% (w/v) PEG 4000 (Sigma) in 50mM Tris- 
HCl, pH 8 and 50mM CaCl 2 ) was added to the tube and 
mixed gently, followed by incubation on ice for lOmin and 
centrifugation at 2000xg for 3min (Megafuge LOR). The 
supernatant was discarded and 1ml of STC buffer (1M 
Sorbitol, 50mM CaCl 2 , 50mM Tris-HCl, pH 8) was added 
to the protoplasts, which were poured onto the surface of a 
prepared PDA plate. The plates were incubated for at 30°C 
for 20hr. 

Batch culture fermentation 

In order to examine the efficiency of knockdown, 
fermentation was carried out with both wild type and 
transformed R. oryzae. Eight batch cultures were set up in 
500ml Erlenmeyer flasks containing 150ml of semi- 
synthetic growth media: 30gm/l glucose, 5gm/l yeast 
extract, 20gm/l (NH 4 ) 2 S0 4 , 20gm/l KH 2 P0 4 , 5gm/l 
MgS0 4 .7H 2 0, 5gm/l CaCl 2 .2H 2 0, 10ml/l trace metal 
solution (3gm/l EDTA (Ci 0 H 14 N 2 Na 2 O 8 -2H 2 O), 0.9gm/l 
CaCl 2 -2H z O, 0.9gm/l ZnS0 4 -7H 2 0, 0.6 gm/1 
FeS0 4 -7H 2 0, 0.2gm/l H3BO3, 0.19gm/l MnCl 2 -4H 2 0, 
0.08gm/l Na 2 Mo0 4 -2H 2 0, 0.06gm/l CoCl 2 -2H 2 0, 
0.06gm/l CuS0 4 -5H 2 0, 0.02gm/l KI), and lml/1 vitamin 
solution (50mg/l d-Biotin, 200mg/l p-amino benzoic acid 
(PABA), lgm/1 Nicotinic acid, lgm/1 Ca-Panthothenate, 
lgm/1 Pyroxidine HC1, lgm/1 Thiamine HC1, 25gm/l m- 
Inositol, pH 5.5). The flasks were inoculated with 20ml/l 
spore solution and cultivated in a shaker incubator at 30°C 
with 140rpm for 62hr. 

HPLC 

The measurement of ethanol, lactic acid, and glucose was 
performed by High-Performance Liquid Chromatography 
(alliance/waters 2695, Waters, USA). The separation 
system consisted of a solvent delivery system equipped 
with an autosampler, a refractive index detector (410 
differential refractometer, Waters), a dual X absorbance 
detector (2487 Millipore, Milford, USA, Waters), and 
computer software (Waters Empower). A hydrogen-based 
ion-exchange column (Aminex HPX- 87H, Bio-Rad, 
USA) was used, and maintained at 65°C with 5mM H 2 S0 4 
eluent at a flow rate of 0.6ml/min. Refractive index (RJ) 
detector and UV absorbance detector at 210nm were used 
in series. Concentrations of ethanol and glucose were 
determined by using the refractive index (RI) detector, 
while the lactic acid was quantified both with RI and UV 
absorbance detectors. 

RESULTS AND DISCUSSION 

Rationale behind efficient siRNA design 

In order to maximize the silencing effect, choosing an 
optimum region of the target gene is critical. It has already 
been reported that among two different genes (IdhA and 



IdhB), coding for NAD + -dependent L-lactate dehydrogen- 
ases (EC 1.1.1 .27) in R. oryzae, IdhA gene plays the key role 
(Skory, 2000). Therefore, IdhA (Accession number 
AB1 11547) was chosen as the prime target for silencing. 
However, as we decided to inhibited both IdhA and IdhB to 
curtail any bypass activity of IdhB, we set out to discover 
regions of extensive homology between the two genes using 
the method of Zhang and Colleagues (Zhang et al, 2000). 
IdhA and IdhB exhibit more than 89% overall nucleotide 
sequence identity and contain no introns. The key sequence 
variation between the two genes lies at position 922-960 in 
IdhB, which IdhA lacks. Using these set of data, the siRNA 
was designed to target nucleotide region 198-222 in the 
IdhA gene, which targeted both IdhA and IdhB mRNAs. 

Efficient R. oryzae protoplast formation 

The cell wall of R. oryzae is mainly composed of chitin 
fibers cemented by chitosan and protein or peptides 
scattered in a mosaic manner (Zamani et al, 2008). 
Therefore, purified lytic enzymes, protease chitosanase, 
and chitinase were necessary for digestion of the cell wall. 
In this study, two different enzymes were investigated for 
degrading the cell wall of the R. oryzae; one based on 
Lysozyme from chicken egg white (Sigma-Aldrich, Fluka 
Biochemika, Switzerland) and another one based on 
Yatalase from culture supernatants of Corynebacterium 
(Takara Bio, Japan). The protoplast was then prepared, and 
the amount of protoplast and formation of mycelia were 
examined accordingly. The results revealed that Lysozyme 
from chicken egg white (Sigma-Aldrich, Fluka 
Biochemika, Switzerland) yielded less protoplast, 
however, growth of R. oryzae on the agar plate was faster 
in those lysed by Lysozyme (data not shown). Formation of 
R. oryzae protoplast by Lysozyme suggests that Lysozyme 
can break down the bond between N-acetyl-D-glucosamine 
residues in chitin present in the cell wall of R. oryzae. 

Effect of IdhA gene silencing on lactic acid and ethanol 
production 

After protoplast preparation by using Lysozyme (lOmg/ml 
in digestion buffer), R. oryzae transformation with siRNA 
was carried out using the calcium chloride/polyethylene 
glycol (CaCl 2 /PEG) method. Different transformation 
systems, such as, electroporation, lithium acetate method, 
microparticle bombardment, and Agrobacterium 
tumefaciens-mediated transformation (ATMT) 
(Nakayashiki, 2008), have been tested in a variety of 
fungi. However in this experiment, the CaCl 2 /PEG 
transformation method was used with minor modifications 
in mucorales fungus R. oryzae. This method was examined 
previously for ascomycetes fungus Neurospora crassa 
(Vollmer and Yanofsky, 1986), Magnaporthe oryzae, and 
Colletotrichum lagenarium (Nakayashiki et al, 2005). The 
results of knockdown were investigated by running eight 
batch cultivations. The product yields obtained for all 
experiments are summarized in Table 1 . 

In six siRNA- inhibited samples, lactic acid production was 
reduced in comparison with the wild type strain. The 
average yield of lactic acid production was O.Olgm/gm in 
inhibited samples, while lactic acid production by wild type 
R. oryzae was 0.07gm/gm. It represents 85.7% (gm/gm) 
reduction in lactic acid production (Table 2). These results 
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Table 1. Product yields from anaerobic batch fermentation in Knockdown and wild type samples of R. oryzae 



Samples 


Lactic acid 
(mg/gm) 


Ethanol 

(mg/gm) 


Glycerol 
(mg/gm) 


Succinic acid 
(mg/gm) 


Pyruvic acid 
(mg/gm) 


KrwM'krinwn U\ 

1U1UI.IVUU TT 11 TT I 


7.7 


450 


72 


0.2 


3.3 


Krwu'lfHftwn HI 

1\IIU *~ IH1 \r TT U TT« 


9.8 


447 


72 


0.6 


2.9 


Knockdown $3 


6.5 


452 


73 


0.2 


3.1 


Knockdown #4 


28.0 


446 


69 


0.2 


2.5 


Knockdown #5 


13.5 


437 


70 


0.1 


2.7 


Knockdown #6 


6.8 


442 


73 


0.2 


2.9 


Wild type#l 


70.9 


411 


64 


0.0 


1.9 


Wild type#2 


66.7 


377 


57 


0.0 


1.7 



Table 2. Average product yields from anaerobic batch cultivations of R. oryzae 



Samples 


Initial 
glucose 
(gm/1) 


Lactic acid 
(mg/gm) 


Ethanol 

(mg/gm) 


Glycerol 

(mg/gm) 


Succinic 

acid 
(mg/gm) 


Pyruvic 

acid 
(mg/gm) 


Biomass 
(mg/gm) 


Carbon 
Balance 


Standard 
Deviation 


siRNA 
Inhibited 


30 


12.05 


445.67 


71.50 


0.25 


2.90 


47.83 


1.057 


0.008 


Wild type 


30 


68.80 


394.00 


60.50 


0.00 


1.80 


53.50 


1.041 


0.038 



indicate noticeable inhibition of lactic acid genes. The 
HPLC results also show that ethanol production increased 
after siRNA-mediated silencing of the Idh genes. The 
average ethanol production increased from 11.82gm/l in 
wild types to 13.37gm/l in knockdown samples (Table 1, 
Figure 1). It corresponds to an average increase in the 
yield of ethanol from glucose from 0.39gm/gm in wild 
types to 0.45gm/gm in the siRNA-inhibited samples 
(15.4% (gm/gm)). 

Effect of RNA silencing on other metabolites of glucose 
metabolism pathway 

Lactate dehydrogenase (LDH) catalyzes conversion of 
pyruvate to lactate. Consequently, the silencing of IdhA, 
i.e., the corresponding gene in lactic acid production, can 
affect the pathway of glucose metabolism in R. oryzae 
(Figure 2). Suppression of lactic acid production in this 
pathway triggers pyruvate to mainly go through the 
production of other metabolites. Moreover, after RNA 
silencing of the lactic acid and ethanol, the yields of 
glycerol, succinic acid, pyruvate, and biomass were 
affected too (Table 1). The yields of glycerol, succinic 
acid, and pyruvate increased after the silencing of IdhA 
gene. The biomass concentration was decreased after the 
silencing of the IdhA gene. The product yields obtained for 
anaerobic batch cultivation of R. oiyzae with initial 
concentration of 30gm/l of glucose as a sole carbon source 
are shown in Table 2. 

Wilcoxon rank sum test 

In order to investigate whether the silencing of the IdhA 
gene resulted in significant differences between the 



production of lactic acid and ethanol in wild type and 
knockdown, respectively, the samples were tested using 
the R statistical program (Version 2.11.1) for one-sided 
Wilcoxon rank sum test (Mann- Whitney U-test), for two 
independent samples as they could not be considered to be 
normally distributed. Wilcoxon rank sum test showed that 
the silenced IdhA gene in R. oryzae could with statistical 
significance, on the 5% level, both reduce the lactic acid 
production (p-value=0.033) and increase the ethanol 
production (p-value=0.036) in knockdowns compared with 
the wild types. 

It is evident that chemically-synthesised siRNAs cannot 
be used on industrial scale to produce ethanol and other 
metabolites due to the prohibitive costs of their 
synthesis. Therefore, the aim of this study was to 
demonstrate that gene silencing by RNAi is indeed 
effective in R. oryzae and can lead to improved ethanol 
yields. We now aim to develop genetically-modified R. 
oryzae strains expressing endogenous shRNAs targeting 
Idh genes for industrial applications. We anticipate that 
the data obtained this study will be of considerable 
value to this future work. 
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Figure 1. Obtained results from HPLC. A. Lactic acid production in Knockdown and wild types Rhizopus oryzae detected by UV in 
HPLC. B. Ethanol production in Knockdown and wild types Rhizopus oryzae detected by RI in HPLC. (WT: Wild type; #1-#1: 
siRNA-treated sample) 
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Figure 2. Glucose metabolism pathway in R. oryzae (modified from Taherzadeh et al, 1996). 
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